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Abstract 

We experimentally demonstrate the nonclassical photon number correlations expected in tripartite con- 
tinuous variable states obtained by parametric processes. Our scheme involves a single nonlinear crystal, 
where two interlinked parametric interactions take place simultaneously, and represents a bright and com- 
pact source of a sub-shot-noise tripartite light field. We analyze the effects of the pump intensities on the 
numbers of detected photons and on the amount of noise reduction in some details, thus demonstrating a 
good agreement between the experimental data and a single-mode theoretical description. 
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I. INTRODUCTION 



Multimode light beams endowed with nonclassical correlations, as those exhibited by multi- 
partite entangled states, represent a resource for quantum technology. They are at the heart of 
enhanced quantum imaging, either ghost imaging or ghost diffraction , and represent a build- 
ing block for the development of an integrated quantum network. In turn, nonlinear interactions 
involving multimode beams of radiation have attracted much attention in the recent years, ei- 
ther to realize all-optical information processing [3] or to generate nonclassical states of light Q4£| . 
Several experimental schemes to generate multimode entangled states have been suggested and 
demonstrated. The first example is provided by the original continuous variable (CV) teleporta- 
tion experiments in Ref. [|5D, where one mode of a twin beam was mixed with a coherent state, 
although no specific analysis was made on the entanglement properties besides the verification 
of teleportation. A similar scheme, where one mode of a twin beam is mixed with the vacuum, 
has been demonstrated and applied to controlled dense coding [Q]. Moreover, a fully inseparable 
three-mode entangled state has been generated and verified by mixing three independent squeezed 
vacuum states in a network of beam splitters |7|]. Recently we suggested and demonstrated a 
compact scheme to realize three-mode entanglement by means of two interlinked x interac- 
tions occurring in a single nonlinear crystal in a type-I non-collinear phase-matching geometry 
IS, 2I. Other schemes involving cascaded interactions have been also analyzed either in period- 



12, 



130. Notice, however, that 



ically poled crystals [10] or in second-order nonlinear ones UllL 
the use of a single nonlinear medium makes the system more compact and robust compared to the 
other schemes that have been suggested and demonstrated so far, in which additional parametric 
sources and linear devices, such as beam splitters, introduce unavoidable losses. Finally, para- 
metric oscillators have been suggested as a source of tripartite signal-idler-pump entanglement in 
triply resonant cavities II 1411 . 

In this paper we experimentally demonstrate the nonclassical photon correlations exhibited by 
tripartite states generated by a single nonlinear crystal, where two interlinked parametric interac- 
tions take place simultaneously. Our scheme realizes a bright and compact source of sub-shot- 
noise three-mode light beams and allows the implementation of simultaneous ghost imaging and 
ghost diffraction protocols with enhanced sensitivity. 

The paper is structured as follows: in the next Section we provide a theoretical description of 
our system and evaluate correlations and noise reduction as a function of the coupling parameters. 
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In Section Un] we describe our experimental apparatus, illustrate the results with focus on nonclas- 
sical photon-number correlations, and analyze the sources of noise in some details. Section [IV] 
closes the paper with some remarks. 



II. THEORETICAL DESCRIPTION 



In our scheme two interlinked interactions, namely a spontaneous parametric downconversion 
process and a sum-frequency generation, take place simultaneously in a single nonlinear crystal. 
In principle, five modes aj are involved in the interactions, two of which, say a 4 and a 5 , are 
non-evolving undepleted pumps and thus are included in the coupling coefficients (parametric 
approximation). The effective Hamiltonian describing the interaction is thus given by 



(1) 



where gi and g 2 are coupling coefficients linearly dependent on the pump fields a 4 and a 5 , respec- 
tively. The earliest studies on the dynamics and the quantum properties of the states realized via 



this Hamiltonian can be traced back to the works in Refs. 115 



urn proi 



16H . The relevance of studying the 



dynamics generated by the above Hamiltonian in details lies in the fact that H int can be realized 
in a variety of different contexts, from quantum optics [floL 

nr 



12 



L20, 



19|] to condensate physics 



210 • The coupling between two optical modes and one vibrational mode of a macroscopic 



object, such as a mirror, has been considered Q22Q and also ions trapped in a cavity have been 
demonstrated to realize the Hamiltonian in Eq. (OQ) for a suitable configuration Q23[|. 

The Hamiltonian admits the constant of motion A(t) = N^t) - N 2 (t) - N 3 (t) = A(0). If we 
take the vacuum |0) = |0)i <g> |0) 2 ® |0) 3 as the initial state, we have N t {t) = N 2 (t) + N 3 (t) Vt, 
being Nj(t) = (aj(t)o/(t)) the mean number of photons in the j-th mode. Under these hypotheses 
the evolved state |T) = exp{— iHi nt t}\0) may be written as 



N n*/2 N r/2 



(m + r)\ 



m\r\ 



m + r,m, r) 



(2) 



+ A^ 1 )( 1+m+r )/ 2 

where we omitted the time dependence of Nj. As a matter of fact the state in Eq. ([2]) is a fully 
inseparable three-mode Gaussian state [24], i.e. a state that is inseparable with respect to any 
grouping of the modes, thus permitting realizations of truly tripartite quantum protocols such as 
conditional twin-beam generation and telecloning [Il7l ll8n. The mean numbers of photons that 
appear in Eq. ([2]) can be obtained by the Heisenberg evolution of the field operators. In particular, 



by introducing f2 = ^\g 2 \ 2 ~~ l^ 1 ! 2 we nave Ni = N 2 + N 3 and 

N 2 = \^pl [cosfit-lf Ar 3 = ^sm 2 (fit). (3) 

We see that when \g 2 \ 2 > \gi\ 2 the dynamics is oscillatory; viceversa, when \gi\ 2 > \g 2 \ 2 we find 
an exponential behavior. 

The above description of the system has been derived under the hypothesis of perfect frequency- 
matching and phase-matching conditions among single-mode fields and the time t appearing in 
Eqs. © represents the interaction time inside the crystal. In this case we did not need to take 
into account the existence of temporal modes and spatial coherence areas. On the other hand, 
if the pump fields are pulsed, the generated fields are temporally multimode [25]. Moreover, 
in a non-collinear interaction geometry, the momentum conservation in the transverse direction 
can be satisfied in more than one way. Thus coherence areas exist [24, whose angles of 
divergence depend on several parameters, such as the pumps intensities, the distance from the 
collinear interaction geometry and the wavelengths of the generated fields. It is interesting to 
point out that in the CV regime the demonstration of the entangled nature of the state in Eq. © 
critically depends on the correct collection of these coherence areas [27]. In fact, collecting light 
from more than a single coherence area corresponds to the introduction of spurious light, while 
collecting less than a coherence area determines a loss of information, which is detrimental to the 
investigations of the nonclassical properties. In addition, we have to select a triplet of areas as 
there is a one-to-one correspondence between the coherence areas in each field. To achieve such a 
selection we can apply a criterion which represents a necessary but not sufficient condition, based 
on the study of the correlation in the number of photons. In fact, due to the constant of motion, 
the state in Eq. is endowed with perfect correlations in the number of photons. The three-mode 
photon-number distribution is given by 

P.(n mr)-5 ^ (m + r) ' (4) 

r T {n, m, r) - d n>m+r ^ - N ^ 1+rn+r mW , W 

from which we can derive the photon-number correlation coefficients between the components 
of the entangled state. In particular, due to the conservation law, we expect the existence of 
strong intensity correlations between the number of photon n% and the sum of the other two, say 
n 2 + n 3 . In order to quantify correlations we denote by ^(nj, rik) = (njUk) — ( n j)( n k) and 
<j2 ( n i) = ( n j 2 ) ~ ( n j) 2 me covariance and the variance of the number of photons, respectively, 



4 



and introduce the correlation coefficients as follows 

7K-,n fc ) 
a(nj)a(n k ) 

Upon exploiting Eq. © we have that the correlation coefficient £1,2+3 is identically equal to one, 
independently of the number of photons generated by the interlinked interactions. On the other 
hand, for the partial photon-number correlations we obtain expressions that do depend on the mean 
number of photons involved. Upon writing N k = PkN where (3\ = (3 2 + P3 and N is the total 
number of photons of the state we have 



N k (l + N{) n»i _ fa-fa 
' J ' A * Xl NA1 + N,) ~ 2AWV W 



, N 2 N 3 N»l _ fc+fc 

( ' 2 "' W (1 + N 2 )(1 + N 3 ) ~ 2(3 2 (3 3 N 1 } 

where from now on k = 2, 3. As the detectors we used to perform the correlation measurements 
are not ideal, we have to rewrite the expressions of the correlation coefficients by taking into 
account the non-unit quantum efficiency of the detection apparatus. The probability operator- 
valued measure (POVM) of each detector, describing the statistics of detected photons, is given by 
a Bernoullian convolution of the ideal number operator spectral measure 



m 



.1 



with j = 1, 2, 3. Equation © can be exploited to calculate the expressions of mean number, (rrij), 
and variance, a 2 (rrij), of the detected photons rrij in terms of the mean number of the photons rij 



and of its variance a (rij) Q3JJ] 



Mj = (rrij) = r ]j (n J )=V j N J (9) 
a\m 3 ) = r,y{n 3 ) + Vj (l - Vj )N, 

We notice that, in general, the statistical distribution of the number of detected photons is different 
from that of the number of photons. Nevertheless, the correlation coefficients, e m , calculated 
for the detected photons can also assume high values; in particular, the correlation coefficient 
calculated between mi and the sum m 2 + m 3 reads as follows 

e m _ Va + N,) AT>>1 _ 1_-Vj_ 



where we have assumed that all the detectors have the same quantum efficiency r\. In turn, the 
partial correlations are given by 

and approach unit value for large N values. 

As a matter of fact a large value of the correlation indices is not sufficient to discriminate 
between quantum and classical correlations [31]. A trivial example is given by the mixture g = 
J2 nmr Pr( n , m,r)\n){n\ <S> \m)(m\ <g> \r){r\, which, with Pt(ti, m, r) given as in Eq. (HI), exhibits 
the same correlations of the state |T). A more realistic example is provided by the tripartite 
state generated by sending a thermal state on two subsequent beam-splitters, whose second port is 
unexcited: the state is classical and shows large intensity correlations, approaching unit value for 



large mean photon numbers [28]. 



In order to obtain a proper marker of nonclassicality we may take into account the difference 



photocurrents dj ik = rrij — m k tl% and build the so-called noise reduction factor 

R k _ g2 fe) (13) 

which is smaller than one for nonclassically correlated states. Note also that, for states generated 
by the Hamiltonian in Eq.©, the existence of sub-shot noise photon-number correlations is a 
sufficient condition for entanglement, i.e. the condition of negative partial transpose is subsumed 
by the condition of sub-shot noise correlations [3]. By using Eqs. © we may write 



T) \(y 2 {n 3 ) +a 2 (n k ) - 2 j(nj,n k )] 
(rij) + (n k ) 

for the noise reduction of bipartite correlations whereas, for the difference photocurrent between 
the mode ai and the sum of the other two modes, we have R = i?i 2+3 



77 

R = 1 - v + - 



E P ° 2 (n p ) + 2T({n k }) 



where 



(15) 



r ({™fc}) = 7(^2, n 3 ) - 7(m, n 2 ) - 7(711, n 3 ) . (16) 

For the state in Eq. © we have 

R=l-V, (17) 



which shows that state |T) exhibits nonclassical tripartite correlations for any value of the mean 
number of photons. Besides, Eq. (fTTT) says that the noise reduction can be detected for any value 
of the quantum efficiency 77. The corresponding bipartite quantities read as follows 

and say that the correlations between modes a 2 and a 3 are always classical whereas the correlations 
between mode a x and either mode a 2 or mode a 3 may be nonclassical in certain regimes. More 



specifically, we have R i:k < 1 if N± < N k + \ / 2N k . Since N% = N 2 + N 3 we may have both 
the noise reduction parameters below the classical threshold only for an overall energy of the state 

N 1 + N 2 + N 3 < 4. 



HI. EXPERIMENT 

The experimental scheme used to generate the nonclassical state of Eq. © is depicted in Fig.Q] 
The harmonics of a continuous-wave mode-locked Nd:YLF laser regeneratively amplified at a 
repetition rate of 500 Hz (High Q Laser Production, Hohenems, Austria) provide the two pump 
fields. In particular, the third harmonic pulse at 349 nm (~ 4.45 ps pulse-duration) is exploited 
as the pump field ct 4 in the downconversion process, whereas the fundamental pulse at 1047 nm 
(~ 7.7 ps pulse-duration) is used as the pump field a 5 in the upconversion process. The two 
processes must simultaneously satisfy energy-matching {uj a = tu 1 + cu 3 , uj 2 = uj 3 + cu 5 ) and phase- 
matching (k| = k° + kg, k 2 = kg + kg) conditions, in which Uj are the angular frequencies, k,- 
are the wavevectors and suffixes o, e indicate ordinary and extraordinary field polarizations. As 
depicted in Fig. [2l we set the pump-field a 4 direction so that the wavevector k 4 is normal to the 
crystal entrance face and propagates along the z-axis of the medium. We also align the wavevector 
k 5 of the other pump field a 5 in the plane (y, z) containing the optical axis (OA) of the crystal and 
the wavevector k 4 . The nonlinear medium is a /5-BaB 2 4 crystal (BBO, Fujian Castech Crystals, 
China, 10 mm x 10 mm cross section, 4 mm thickness) cut for type-I interaction ($ cut = 38.4 
deg), into which both pumps are strongly focused. Typical intensity values of the pumps were ~ 5 
GW/cm 2 for a 4 and ~ 2 GW/cm 2 for 05. The required superposition in time of the two pumps is 
obtained by a variable delay line. 

With reference to Fig.[2l we indicate as dj the angles in the plane (y, z) formed by each wavevector 
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with k 4 and as (3j the angles of each wavevector with respect to this plane. For the experimental 
realization of the interaction scheme we choose the solutions in the plane (y, z), thus (3j = for 
j = 1 — 3: in particular, we sent the pump field a 5 at an external angle $5 : ext = —24.47 deg with 
respect to the other pump field a 4 . Under these hypotheses, for Ai = 632.8 nm, A 2 = 446.4 nm 
and A 3 = 778.2 nm, we calculated the following external interaction angles with respect to the 
pump field a 4 : $i, e xt = -9.78 deg, i? 2 ,eart = -3.25 deg and $ 3>ext = +12.06 deg 

The preliminary use of a He:Ne laser as the seed field allowed us to position three pin-holes 
on the path of the three generated fields in such a way that then, when operating the system from 
vacuum (i.e. in the absence of any seed fields), we could collect a triplet of coherence areas. 
Distances and sizes of the pin-holes were chosen by searching for the condition of maximum 
intensity correlations between the generated fields 112 811 . In fact, as shown in Section HH we expect 
strong intensity correlations not only between the number of detected photons mi and the sum of 
the other two, but also between m x and m 2 , m 2 and m 3 and mi and m 3 . By applying this criterion, 
we finally decided to put two pin-holes of 30 fim diameter at distances d\ = 60 cm and d 3 = 49 
cm from the BBO along the path of the signal beam at 632.8 nm and of the idler beam at 778.2 
nm, respectively. The two different distances were chosen to compensate for the difference in the 
divergence of signal and idler due to their wavelengths 112411 . Moreover, as the beam at 446.4 nm 
has a divergence smaller than those of the other two fields, we selected it by means of a 50 /xm 
diameter pin-hole placed at a distance d 2 = 141.5 cm from the crystal. 

The light, suitably filtered by means of bandpass filters located in front of each pin-hole, was 
focused on each detector by a lens (f\ = / 3 = 25 mm, f 2 = 10 mm). Since we performed 
measurements in the macroscopic intensity regime (more than 1000 photons per coherence area), 
we used three p-i-n photodiodes (two, D x 2 in Fig. [T] S5973-02 and one, D 3 , S3883, Hamamatsu, 
Japan) as the detectors. In order to obtain the same overall detection efficiency (bandpass filter 
plus detector) on the three arms, we put two adjustable neutral-density filters in the pathways 
of a 2 and a 3 , thus obtaining the same value r] = 0.28 on the three arms. The current output 
of the detectors was amplified by means of two low-noise charge- sensitive pre- amplifiers (CR- 
110, Cremat, Watertown, MA) followed by two amplifiers (CR-200-4 /is, Cremat). We connected 
the detectors D 2 and D 3 to the same amplifier device by means of a T-adapter. The two amplified 
outputs were then integrated by synchronous gated-integrators (SGI in Fig.Q]) operating in external 
trigger modality (SR250, Stanford Research Systems, Palo Alto, CA). The voltage outputs were 
then sampled, digitized by a 12-bit converter (AT-MIO-16E-1, DAQ National Instruments) and 
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recorded by a computer. 

In the following we discuss the measurements of the intensities of field a\ and of the sum 
a 2 + a 3 as, according to Eqs. (TT7T)-(fT9l), we expect a nonclassical behavior. Partial measurements 
performed by alternatively blocking the light impinging on the detectors D 2 and _D 3 were not 
very reliable as the numbers of detected photons on the two fields separately were too close to 
the electronic noise of the detection chain. This is an important drawback as, for all calculations 
that follow based on experimental data, we must take into account the electronic noise that we 
measured in the absence of light I27L 

As the pump fields are pulsed and their duration is longer than the characteristic time of the 
nonlinear processes ,the distributions of the detected photons collected by the pin-holes are tem- 
porally multimode |25|]. The same is also true for the statistical distribution of the sum m 2 + m 3 . 



Moreover these distributions should be characterized by the same number of modes [24] 



From the experimental point of view, the main difficulty to be overcome was the correct selec- 
tion of a triplet of coherence areas. In fact, in the CV domain we have to avoid spurious light that 
could be detrimental to the experimental results; moreover, the interaction scheme presented here 
involves not only two generated fields, but also a third one, which obviously makes the detection 
more critical. Finally, we have two pump fields instead of one and in particular we are not able 
to exactly measure the effective portions of them that interact into the crystal. In spite of all these 
difficulties, we characterized the state produced by the interlinked interactions and in particular we 
proved its quantum nature by performing sub-shot noise photon-number correlation measurements 
as a function of the pump fields intensities. In fact, as remarked in SectionUIl the evaluation of the 
noise reduction factor R for the distribution of the difference photocurrent d — m\ — (m 2 + m 3 ) 
provides a sufficient condition in order to test the quantum nature of the generated state. 

We firstly investigated the evolution of the mean number of photons as a function of the inten- 
sity of one of the two pumps by keeping fixed the intensity of the other one Q29[|. In fact, if on 
one hand this analysis allows us to verify that the mean number of photons does not depend on the 
correct selection of the coherence areas, on the other one it is essential for the determination of the 
effective values of the pump fields intensities from the fitting curves. 

As a first check we studied the evolution of the mean number of detected photons, Mi and M 2 + 
M 3 , as a function of the intensity of field a 4 for a fixed value of the intensity of field a 5 . Note that 



temporal evolution in Eqs. © is transformed into spatial evolution by identifying a/I^I 2 — \gi\ 2 t 
with A/I72I 2 — |7i| 2 <2> z being the effective interaction length 113011 . In the experimental condition 
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each Mj represents the total mean number of photons detected beyond each pin-hole; actually, it 
can be expressed as Mj = n(m,j), where \i is the number of temporal modes and (mj) the average 
population of each mode. To vary the intensity of field a 4 , we changed its energy by means 
of an adjustable neutral-density filter. For each energy value, measured by means of a movable 
thermal detector (D A in Fig. [TJ mod. 03A-P-CAL-SH, Ophir Optronics Ltd., Jerusalem, Israel), 
we measured the mean number of photons by averaging over 50000 subsequent laser shots. In 
Fig. Oa), we show the measured values of M\ and M 2 + M 3 as functions of |7i| 2 , for a fixed value 
of |7 2 | 2 . Note that |7i| 2 oc E 4 / '(7rr|/za> 4 r 4 ), E 4 being the pulse energy of field a 4 , r 4 the pulse 
duration and r 4 the beam radius. The experimental data are displayed together with their common 
fitting curve, obtained from Eqs. © with |7 2 | 2 as the parameter and | '-y 1 1 2 as the variable. In this 
case we get \ j 2 1 2 = 8.17 x 10 5 m~ 2 and (71 1 2 in the range 1.86 x 10 6 - 2.17 x 10 6 m" 2 . Note 
that, as expected, the experimental data satisfy the photon-number conservation law as they are 
almost superimposed. The best fitting curve has been obtained allowing a slight difference in the 
quantum efficiency values of the detection chains and finding the values from the conservation law. 
We found 771 = 0.31 and i] sum = 0.28. The difference is within the error justified by the tolerance 
of the pin-holes sizes (0i = 3 = 30 ± 2 fxm and 2 = 50 ± 3 fim) and it is also justified by 
possible imperfections in the positioning of the pin-holes at the right distances from the crystal. 

As a second check, we studied the evolution of Mi and M 2 + M 3 as a function of the intensity 
of field a 5 , by keeping the intensity of field a 4 fixed. To change the energy of field a 5 we placed 
a half-wave plate on the pathway of the infrared pump field. A movable thin-film plate polarizer 
was used to measure the energy fraction corresponding to the ordinarily polarized component of 
the field for each step of rotation of the A/ 2 plate. For each energy value, measured by means of 
the thermal detector (_D 5 in Fig. [D, we measured the mean number of photons by averaging over 
50000 subsequent laser shots. In Fig.[3£b), we show the measured values of Mi and M 2 + M 3 as 
functions of |7 2 | 2 , for a fixed value of |7i| 2 . Also in this case, the experimental data are plotted 
together with the fitting curve of the two sets of data obtained from Eqs. ©. Obviously, we 
have to interchange the roles of the pumps: in fact, here |7i| 2 is treated as the parameter and 
I72I 2 as the variable. In particular, we obtained |7i| 2 = 1.52 x 10 6 m~ 2 and |7 2 | 2 in the range 
(1.97 x 10 4 — 1.27 x 10 5 ) m~ 2 . Even in this case, the experimental data satisfy the conservation 
law as they are almost superimposed and the optimization of the quantum efficiencies still gives 
very small corrections: % = 0.283 and r] sum = 0.28. 

By exploiting the values of the pump fields intensities obtained from the fitting curves, we 
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can investigate the behavior of the correlation coefficient e™ 2+3 (see Eq. (flOl) ) and of the noise 
reduction R (see Eq. (fT7T) ). First of all, in Fig. H] we show the intensity correlation coefficient, (a), 
and the noise reduction, (b), as functions of |7x| 2 by keeping fixed the value of |7 2 | 2 . During these 
measurements the collection areas were kept fixed (same pin-holes located at the same distances 
as above). The variation of the correlation coefficient as a function of the pump field intensity 
through |7i| 2 is indeed not so strong, but the noise reduction factor is critically dependent on the 
changes in the intensity value. In fact, very much as in the case of the twin-beam state [27], 
there is an optimum condition at which R is minimum and, correspondingly, the value of the 
correlation coefficient is maximum. Note that |7i| 2 is larger than |7 2 | 2 in the entire range of 
variation. Moreover, we note that increasing the pump intensity, hence |7i| 2 , also increases the 
size of the coherence areas so that they are only partially transmitted by the pin-holes. On the other 
hand, lowering the pump intensity reduces the size of the coherence areas and allows uncorrected 
light to pass the pin-holes. Note that the values of R corresponding to the selection of more than a 
single coherence area remain quite close to the shot-noise limit as the information contained in the 
area is not lost, but only made more noisy. On the contrary, the selection of only a part of the areas 
causes a loss of information that determines a more remarkable increase of R above the shot-noise 
limit (note the axis break). This result represents an indication of the need of a perfect matching 
of the pin-hole areas in order to obtain sub-shot noise correlations. 

Secondly, we investigated the intensity correlation coefficient and the noise reduction as func- 
tions of 1 72 1 2 by keeping fixed the value of |7i| 2 . Also in this case the collection areas were kept 
fixed by using the same pin-holes as before located at the same distances. The intensity regime 
in which these measurements were performed is different from the previous one as the absolute 
values of the two pump fields intensities are smaller than in the other case. However, |7 2 | 2 is again 
smaller than |7i| 2 in all its range of variation. For all these reasons, the variations in the exper- 
imental values of the correlation coefficient and of the noise reduction are smaller (see Fig. [5]). 
Moreover, the minimum value of R is quite near to the lower limit R m %n = 0.72 fixed by the quan- 
tum efficiency. In fact, the use of less intense pumps reduces the quantity of spurious light that 
can be revealed by the detectors; in addition, the fluctuations of the laser source and the possible 
discrepancy of its photon-number distribution with respect to the ideal Poissonian statistics play a 
less important role J3I I. 

As a further investigation, we performed other measurements in order to study how critical 
is the sub-shot noise condition with respect to a slight change in the values of the intensities of 
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the two pump fields. In particular, we verified that it is always possible to choose the pumps in 
such a way that only micro-metric adjustments of the pin-holes positions are necessary to select 
the coherence areas. In Fig. [6] we show a number of sub-shot noise measurements obtained for 
different pairs of pump values. As we can see, not all the measurements reach the optimum 
minimum value, R ~ 0.72, due to residual imperfections in the selection of the coherence areas. 
In particular, as we remarked above, this operation is more critical when the intensity values are 
higher because other noise sources become important 113 ill . However, we want to emphasize that 
there are many possible choices of the pumps values that allow us to perform sub-shot noise 
measurements thus demonstrating that our scheme is particularly versatile and useful for several 
applications in different photon-number regimes. 



IV. CONCLUSIONS AND OUTLOOKS 



In conclusion, we have presented the experimental realization of an entangled state that in- 
volves three modes of radiation in the macroscopic regime. We verified the quantum nature of the 
state produced by our all-optical interaction scheme by means of sub-shot noise photon-number 
correlations, which also subsumes the inseparability condition. In particular, we investigated how 
critical is the sub-shot noise condition by studying its dependence on the intensities of the two 
pump fields. In spite of the difficulties in measuring the light of triplet coherence areas correctly 
and in avoiding the detection of spurious light, we obtained quite relevant results that could be 
further optimized. In the immediate future we plan to use three acquisition chains instead of only 
two to separately and simultaneously detect the three fields. Moreover, in order to reduce the noise 
that is detrimental to the shot-noise reduction factor we still intend to operate in the macroscopic 
regime, but with lower numbers of photons and to use hybrid photodetectors, which are endowed 
with a reasonable quantum efficiency (77 ~ 0.4) and a linear response in the mesoscopic regime 
(up to a few hundreds of detected photons), instead of the p-i-n photodiodes. We also plan to mod- 
ify our collection system by using optical fibers in order to avoid spurious light and to minimize 
the uncertainty in the collection areas. The experimental improvements would make the whole 
system more easily controllable and suitable for several applications, such as the production of 
conditional twin-beam states and the generation of quasi-Fock states with a number of photons 
sensibly greater than one. Overall, our system represents a robust and tunable scheme to obtain 
nonclassical photon number correlations in tripartite CV systems thus allowing the simultaneous 
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realization of ghost-imaging and ghost-diffraction with enhanced sensitivity. 
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FIG. 1: Scheme of the experimental setup: BBO, nonlinear crystal; NF, variable neutral-density filter; A/2, 
half-wave plate; TFP, thin-film plate polarizer; P1-3, pin-holes; /i_5,5', lenses; Di_3, p-i-n photodiodes; 
D45, thermal detectors; M, Aluminum mirrors; PRE+AMP, low-noise charge-sensitive pre-amplifiers fol- 
lowed by amplifiers; SGI, synchronous gated-integrator; ADC+PC, computer integrated digitizer. 




FIG. 2: Scheme of the phase-matched interlinked interactions: (x, y) -plane coincides with the crystal en- 
trance face; a, tuning angle; /3/s, angles to (y, z)-plane; i?/s, angles on the (y, z)-plane; ip, angle to the 
optical axis (OA). 
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FIG. 3: (a) Evolution of the mean numbers of detected photons as a function of 1 71 1 2 which is proportional 
to the intensity of field 04 for I72I 2 = 8.17 x 10 5 m~ 2 . Black circles: measured values of Mi; grey 
triangles: measured values of M2 + M3; solid straight line: fitting curve, (b) Evolution of the mean 
numbers of detected photons as a function of I72I 2 which is proportional to the intensity of field 05 for 
|7i| 2 = 1.52 x 10 6 m~ 2 . Black circles: measured values of Mi; grey triangles: measured values of 
M2 + M3; solid straight line: fitting curve. 
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FIG. 4: (a) Intensity correlation coefficient and (b) quantum noise reduction R (note the axis break) as 
functions of | 7 i| 2 for | 72 | 2 = 8.17 x 10 5 m" 2 . 
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FIG. 5: (a) Intensity correlation coefficient and (b) quantum noise reduction R (note the axis break) as 
functions of | 72 | 2 for | 7 i| 2 = 1.52 x 10 6 m' 2 . 
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